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Global Sea Level Trend in the Past Century 


Abstract. Data derived from tide-gauge stations throughout the world indicate 
that the mean sea level rose by about 12 centimeters in the past century. The sea 
level change has a high correlation with the trend of global surface air temperature. 
A large part of the sea level rise can be accounted for in terms of the thermal 
expansion of the upper layers of the ocean. The results also represent weak indirect 
evidence for a net melting of the continental ice sheets. 


Sea level change is of current interest 
because of its possible sensitivity to cli- 
mate change. It has been suggested, for 
example, that global warming due to 
increasing atmospheric CO, could melt 
the marine West Antarctic ice sheet, 
raising the global sea level 5 to6m(/). A 
sea level rise of as little as 15 cm may 
double the probability of damaging storm 
surges on the coast of Britain (2). Such a 
rise would also cause substantial beach 
erosion and the intrusion of seawater 
into low-lying areas that are now fresh- 
water regions. 

Many processes affect the sea level 
position measured on shorelines. Among 
the most important are eustatic sea level 
changes due to changes in the ocean 
water volume, caused mainly by the 
melting or growth of ice sheets, and 
isostatic adjustments of the earth’s crust, 
caused mainly by ice sheet growth or 
decay and the associated change in the 
ocean water mass (3, 4). Tectonic move- 
ment and river sedimentation can gener- 
ate local sea level trends comparable to 
eustatic and isostatic changes, as can 
changes in ocean currents and prevailing 
winds, although such trends are of limit- 
ed duration. 

We used tide-gauge measurements to 
estimate global sea level change in the 
past century. Data from more than 700 
stations were obtained from the Institute 
for Oceanographic Science, Birkenhead, 
England. We excluded stations with rec- 
ords shorter that 20 years (a majority of 
the stations) and stations in seismically 
active areas such as the Pacific coast of 
Japan and in rapidly subsiding localities 
such as Galveston, Texas, and the Mis- 
sissippi delta. The remaining 193 stations 
were divided into 14 regions on the basis 
of geographic proximity and the expect- 
ed similarity of isostatic or tectonic be- 
havior. 

We reduced the individual station rec- 
ords to a common reference point by 
fitting a least-squares regression line to 
sea level as a function of time and by 
defining the zero point to be the value of 
the regression curve for 1940. The annu- 
al mean sea level curves for stations 
within a geographical region were then 
averaged to yield a mean sea level curve 
for each region. We obtained the global 
mean sea level curve by averaging the 


regional mean sea level curves, weight- 
ing each region equally (excluding the 
isostatically uplifting region of Scandina- 
via). Sea level trends, obtained by fitting 
a regression line through the mean sea 
level curves, are summarized in Table 1. 

We also attempted to remove the long- 
term (usually 6000-year) sea level trends 
from the station data in order to obtain 
short-term sea level fluctuations, which 
are perhaps more appropriate for corre- 
lation with global climate variations in 
the past century. The cause of the long- 
term trend is uncertain. It has been ar- 
gued that as much as 90 percent of it is 
residual isostatic uplift of continents due 
to the removal of the Wisconsin ice 
sheets (4). However, the long-term trend 
may contain a eustatic component, for 
example, due to a change in volume of 
the Antarctic or Greenland ice sheets. 

The estimates for long-term sea level 
change are based on “C dating of mea- 
sured positions of shoreline indicators in 
the geologic records, for example, mol- 
lusks, corals, and brackish-water peats 
(5). The 6000-year time interval was cho- 
sen to be as large as possible without 
approaching the period of the North 
American and European ice sheets, thus 
minimizing the effect of errors in the 
estimated sea level trend. However, for 
Scandinavia, where there has been a 
high rate of isostatic uplift since the last 
deglaciation, we used as a time interval 
the last 2500 years to improve the likeli- 
hood of obtaining a linear trend applica- 
ble to the past century. We corrected the 
dating of sea level positions to use the 
recent 5730-year estimate for C half- 
life (6), rather than the 5570-year stan- 
dard that had been used for most of the 
records (5), and to account for past at- 
mospheric '4C fluctuations (7). 

The raw global sea level trend that we 
obtain for the past century is a rise of 
about 12 cm per century. After subtrac- 
tion of the long-term trend, the result is a 
rise of 10 cm per century or 1 mm/year. 
The inferred long-term trend of 2 cm per 
century is very small as compared to the 
long-term trends of 1 m per century or 
more that are common at times of conti- 
nental ice sheet growth or decay (3). 
Evidently the past few millennia have 
remained too warm to permit ice sheet 
formation on the North American or 
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Eurasian continents but too cold for sub- 
stantial melting of the Greenland or Ant- 
arctic ice sheets. 

We find that sea level rose in the past 
century in every geographical region ex- 
cept Scandinavia, and, after correction 
for long-term trends, sea level rose in 
every region except the west coast of 
South America where the change is 
smaller than the uncertainty based on the 
95 percent confidence limit. We thus 
believe that this sea level rise is a true 
global trend and not, for example, a 
result of some regional variation in the 
geoid. 

The sea level trend we find is similar to 
that obtained by Gutenberg [1.1 mm/year 
for the period 1807 to 1939 (8)], Fair- 
bridge and Krebs [1.2 mm/year for 1900 
to 1950 (9)], and Lisitzin [1.12 mm/year 
for 1807 to- 1943 (/0)], even though their 
studies were based on a much smaller 
number of tide-gauge stations. However, 
this rise is much less than the value 
recently reported by Emery [3 mm/year 
for 1935 to 1975 (/1)]. Emery included 
stations in regions of known strong local 
uplift and subsidence (for example, 
Scandinavia and the east coast of Japan) 
and weighted each station equally; the 
result was that Scandinavia, Japan, and 
the east coast of the United States were 
heavily weighted. Emery (//) also re- 
ported a large increase in sea level (7 
mm/year) in the decade 1966 to 1975 on 
the southeast coast of North America, 
and Wanless and Harlem (/2) claim that 
sea level rose 10 to 14 cm in the past 
decade in south Florida; their result is 
based on the displacement of intertidal 
organisms on seawalls and bridge pil- 
ings. Our data (Fig. 1) show the steep 
rise in that region for 1966 to 1973, but 
the rise was not global and sea level 
decreased in the eastern United States 
for the period 1973 to 1977. Evidently 
long-term trends cannot be estimated 
from changes over 5- to 10-year periods 
for a single region. Even the global aver- 
age curve (Fig. 2) has notable short-term 
variability. 

The global sea level trend for the past 
century has some similarity to the trend 
in global surface air temperature (/3). To 
quantify this, we computed the correla- 
tion coefficient between our global sea 
level curve and the global temperature 
curve of Hansen et al. (13), obtaining 0.6 
when the annual mean curves are used 
for both quantities and 0.8 when the 5- 
year running means are used. Most of 
the positive correlation arises from the 
general increase in both sea level and 
temperature. Since one might expect a 
time lag between temperature change 
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and sea level, we fitted a linear relation 
between our sea level curve and the 
global temperature trend, 


S) = aT — to) + b (1) 


where S and T are the 5-year means of 
global sea level and temperature, respec- 
tively, and ż is time. The parameters a 
and b were obtained by least-squares 
linear regression, and the time lag tọ was 
chosen to minimize the variance be- 
tween Eq. | and the sea level curve. The 
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Fig. 1. Regional mean sea level trends. The 
heavy lines are 5-year running means. Long- 
range (6000-year) trends have been subtract- 
ed. 


results were a = 16 cm K™!, b = 0.3 
cm, and tọ = 18 years. 

The sea level and temperature records 
are too short to allow much significance 
to be attached to this relationship. Nev- 
ertheless, it is interesting that the time 
lag of 18 years is of the order of the 
thermal relaxation time for the upper 
layers of the ocean, that is, the layers 
that are mixed in a time less than or 
comparable to the thermal relaxation 
time. This result suggests that part of the 
sea level rise may be attributable to 
thermal expansion, a possibility we can 
test by using the heat fluxes from the 
one-dimensional model of Hansen et al. 
(13). The two primary parameters or 
“tuning knobs” in that model are the 
equilibrium sensitivity (say, ATeq for 
doubled atmospheric CO.) and the rate 
of mixing of heat into the ocean beneath 
the mixed layer (specified by a diffusion 
coefficient k). Figure 2 illustrates the sea 
level change obtained for three values of 
the model’s equilibrium sensitivity, k be- 
ing constrained in each case to the value 
providing the best agreement with the 
observed global temperature trend of the 
past century. With the commonly ac- 
cepted value for AT, of ~ 3°C, about 
half the observed sea level change is 
accounted for in terms of the thermal 
expansion of seawater. We tried other 
models for the mixing of heat into the 
ocean and obtained similar results. 

We conclude that a large part, but 
probably not all, of the sea level rise of 
the past century is due to thermal expan- 
sion of the upper ocean. These results 
therefore also provide weak evidence for 
a decrease in the volume of the nonocean 
reservoirs of water. Some ground-water 
levels are known to have receded recent- 
ly. For example, the huge Ogalalla reser- 
voir in the high plains of the central 
United States may have dropped by sev- 
eral tens of meters (74), equivalent to a 
drop of a few millimeters of sea level. It 
seems possible that a sea level rise of a 
few centimeters could be accounted for 
in terms of a lowering of global ground- 
water levels. On the other hand, the 
trapping of water behind dams in the past 
century may have reduced sea level by 1 
to 2 cm (/5). Thus the net change in the 
water reservoirs of land areas is proba- 
bly not a major cause of sea level 
change. 

The magnitude of the global sea level 
rise therefore suggests that there has 
been a small net melting of the ice 
sheets. Although the present evidence is 
weak, continuation and refinement of 
this type of analysis is potentially impor- 
tant. As yet, direct observations of ice 
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Table 1. Sea level trends, 1880 to 1980, including correction for long-term (6000-year) trends. 


Sea level trend, 1880 to 1980 


Corrected sea level trend, 1880 to 1980 


Num- Ener 95 percent Num- Linear 95 percent 
Region ber : confidence ber confidence 
trerid phe trend tee 
of (cm/100 limit of (cm/100 limit 
sta- (cm/100 sta- : (cm/100 
tions years) years) tions years) years) 

West coast, North America 16 10 2 1 8 3 

Gulf Coast and Caribbean 6 23 4 4 16 5 

East coast, North America 32 30 2 30 15 2 

Bermuda 1 26 16 1 20 16 

West coast, South America 8 19 31 2 -3 3 

East coast; South America 5 4 11 2 16 11 
Africa 2 32 31 0 

Southern Europe 15 32 2 7 7 2 

West central Europe T 13 2 5 4 2 

Southern Baltic 21 4 2 14 5 2 

Scandinavia 47 —37* 3* 10 10 3 

Asia 9 4 3 2 22 4 
Australia 9 13 3 0 

Pacific Ocean 15 19 3 6 6 4 

Global mean 193 12 1 86 10 1 


*Not included in the global average. 


sheet mass balance are not sufficiently 
accurate to establish even the sign of any 
trend (16). 

A key application of the global sea 
level trend concerns the potential de- 
struction of the marine West Antarctic 
ice sheet (/). It can be argued that the 
ice sheet is not close to disintegration, 
because it survived the Altithermal 
(~ 5000 years ago) when the global 
mean temperature was perhaps 1°C 
warmer than today. However, sea level, 
as well as temperature, must affect the 
ice sheet’s stability. As indicated above; 
sea level has been flat, perhaps even 
rising slowly, over the past 5000 years. 
With the 10-cm rise of the past century, 
sea level must now be at or near its 
highest level since the previous inter- 
glacial, the Eemian 120,000 years ago 
(17). 

Thus it is not inconceivable that the 
situation is near a point at which contin- 
ued warming and rise of sea level could 
cause rapid, highly nonlinear disintegra- 
tion of the ice sheet (J). We should 
emphasize that we have no evidence for 
such a process. Indéed, the sea level 
change we have deduced appears to be 
linear with temperature, and largely a 
result of the thermal expansion of seawa- 
ter. Nevertheless, since sea level is at a 
high point and rising, the West Antarctic 
ice sheet warrants close attention. 

Continued rise of sea level is likely in 
the near future, if predictions of global 
warming (13) are correct. The thermal 
expansion of seawater may raise sea 
level about 20 to 30 cm in the next 70 
years (8); if slow ice sheet melting in- 
creases this by the same factor as in the 
past 100 years, a sea level rise of about 
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40 to 60 cm would occur by 2050. Thus 
we believe that substantial sea level 
change may occur even without rapid 
collapse of the West Antarctic ice sheet. 

There is a clear need for improved 
observations. A direct measure of ice 
sheet growth or decay could be obtained 
from observations of ice sheet areal cov- 
erage and the altitude of the upper sur- 
face of the ice sheets that can be moni- 
tored by satellite. It is also desirable to 
measure the changes in the temperature 
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profile in the ocean; these measurements 
would provide more explicit information 
on heat penetration and would serve as a 
crucial test of ocean-atmosphere climate 
models. 

Finally, there is a need for observa- 
tions and studies of sea level in the 
1980’s. A sharp global warming trend has 
been under way since the mid-1960’s 
(13), and the current growth of atmo- 
spheric CO, and trace gases virtually 
assures that this trend will continue (/9). 


e.e... ATog=5.6°C (k= 2.2 cm? sec!) 


ATeq=2.8°C (k=1.2 cm? sec”) 


— = A Teq=1.4°C (k=0 cm? sec™') 
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Fig. 2. Global mean sea level trend from tide-gauge data and comparison to the thermal 
expansion of the upper ocean obtained from the model of Hansen et al. (13) (see their equation 9 
for the heat flux into the ocean). The radiative forcing used was CO, + volcanoes + sun [figure 
5 in (73)}, but a similar result would be obtained for other forcings that fit the observed global 
temperature trend; AT, is the equilibrium sensitivity of the model for doubled CO,, and & is the 


diffusion coefficient beneath the mixed layer. 
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The sea level response to this warming 
should be carefully determined to aid our 
understanding of the processes and to 
allow early detection of any nonlinear 
response. 
V. GORNITZ 
S. LEBEDEFF 
J. HANSEN 
Goddard Space Flight Center, 
Institute for Space Studies, 
New York 10025 
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represents a source of fast oxygen atoms in 


Venus’ exosphere, and subsequent collisions of oxygen atoms with hydrogen atoms 
lead to escape of about 10’ hydrogen atoms per square centimeter per second. 
Escape of deuterium atoms is negligible, and the ratio of deuterium to hydrogen 
should increase with time. It is suggested that the mass-2 ion observed by Pioneer 
Venus is D*, which implies a ratio of deuterium to hydrogen in the contemporary 
atmosphere of about 1077, an initial ratio of 5 x 107°, and an original H,O 
abundance not less than 800 grams per square centimeter. 


Venus contains quantities of carbon 
and nitrogen similar to Earth, but hydro- 
gen is deficient (7). The abundance of 
water on Venus is about 4.2 g cm’, 
which may be compared to the terrestrial 
value of 2.7 x 10° g cm~*. Walker et al. 
(2) argued that the amount of H,O on 
Venus should have been initially similar 
to that on Earth; they proposed that 
escape of H from a hot steamy atmo- 
sphere played a major role in Venus’ 
early evolution. Lewis (3) thought that 
the H,O content was low from the out- 
set, because Venus formed in a warmer 
region of the solar nebula. 

Analysis of processes influencing the 
current budget of Venus’ hydrogen 
should shed some light on this issue. The 
upper atmosphere contains significant 
quantities of H, measured first by Mari- 
ner 5 (4). The density of H in outer 
regions of a planetary atmosphere pro- 
vides direct information on the exchange 
of hydrogen with the interplanetary me- 
dium, and analysis of data for Earth (5) 
indicates an escape flux for H of 
2.7 x 108 cm~? sec™!. The distribution 
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of H is more complex for Venus, with at 
least two components (6), and interpreta- 
tion is correspondingly ambiguous (7). 

One component of Venus’ exospheric 
H has a scale height of about 300 km, 
consistent with temperatures observed 
by Pioneer Venus (8). The second is 
more extensive, with a scale height of 
about 1000 km. Escape of atoms in the 
cold component is trivial, about 10 cm~? 
sec™!. Escape from the extended distri- 
bution may proceed more readily. The 
escape rate would be 1.5 x 10° cm~? 
sec”! if the velocity distribution were 
characterized by an effective tempera- 
ture of 1000 K as suggested by the ob- 
served scale height. 

Fast atoms comprising the extended 
component must be produced in the exo- 
sphere, either by acceleration of ambient 
H or by exothermic reactions involving 
hydrogen-bearing gases such as H3. Pro- 
posed chemical sources (9, /0) include 
charge transfer of H with solar wind and 
a variety of reactions involving iono- 
spheric ions, for example, reactions 1 
through 6 in Table 1. We propose here an 
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additional source, collisions of fast O 
with thermal H, reaction 7. Energetic O 
atoms, O*, are formed by the recombina- 
tion of exospheric O,* and CO)", reac- 
tions 8 and 9. Recombination of O,* 
proceeds mainly (//) by reaction 10, 
which represents a source for O* with an 
initial speed of 5.6 km sec™!. Oxygen 
atoms formed in reaction 10 can escape 
directly from Mars (12). They are gravi- 
tationafly bound to Venus but can induce 
significant escape of H by momentum 
transfer in reaction 7. 

Consider an elastic collision of O* 
with stationary H and assume that the 
scattering is isotropic. The fraction of 
collisions leading to production of H 
atoms with speeds in excess of the es- 
cape velocity, 10.2 km sec”!, is 6 per- 
cent for a O* velocity of 5.6 km sec”! 
(13). The fraction is increased to approx- 
imately 15 percent if we account for the 
thermal velocity of H in an exosphere of 
temperature 300 K as implied by the cold 
H component observed by Mariner 5 
(6). 

The rate for escape of H due to reac- 
tion 7 on the dayside may be estimated 
as follows. Photoionization above the 
exobase leads to the production of O,* 
and subsequently O*. The rate for pro- 
duction of O* is given by 2/N, where J 
denotes the ionization rate (per second) 
and N is the column density of gas above 
the exobase (per square centimeter). 
Typical values for J range from 5 x 1077 
to about 1.2 x 107°, with higher values 
appropriate for solar conditions during 
the Venus flyby of Mariner 5 (/4). The 
value of N is about 3 x 10! (15), and the 
resulting source for O* is 8 x 108 em~? 
sec”!. The fraction of exospheric O* that 
collides with H depends on the relative 
abundance of H and is approximately 2 
percent for conditions during the Mari- 
ner 5 flight (/6). The corresponding rate 
for escape of H from the dayside is about 
10° cm? sec™!. 

The source of exospheric O* and H* 
may be even larger at night. The night- 
side ionosphere is variable (/7), with 
densities as high as 10° electrons per 
cubic centimeter at altitudes between 
140 and 150 km (/8). The dominant posi- 
tive ion at these altitudes is O,*, pro- 
duced by reaction 11. The Pioneer data 
imply a recombination rate for the ions 
of about 3 x 108 cm? sec! (18, 19), 
which may be supplied either by trans- 
port from the dayside thermosphere or 
by in situ production (7-20). The corre- 
sponding rate for production of nightside 
O* is 6 x 108 cm~? sec |. In contrast to 
the dayside, most of this source is locat- 
ed in the exosphere. The density of H in 
the lower nightside exosphere is 100 
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